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HYPERSONIC  AXISYMMETRIC  WAKES 
INCLUDING  EFFECTS  OF  RATE  CHEMISTRY 

by  M.H.  Bloom  and  M.H.  Steiger 

I  INTRODUCTION 

A  wake  is  the  flow  field  downstream  of  a  body,  and  is  characterized 
by  disturbed  properties  which  eventually  return  to  ambient  atmospheric 
conditions.  Wake  characteristics,  such  as  velocity  defects,  may  be  gene¬ 
rated  by  the  boundary  layer  shed  by  the  body  under  subsonic  or  supersonic 
conditions,  or  may  be  generated  by  the  bow  shock  under  supersonic  condi¬ 
tions,  or  by  the  combined  effects  of  boundary  layer  and  shock.  Evidently 
the  relative  importance  of  the  boundary  layer  and  shock,  insofar  as  wake 
properties  are  concerned,  depends  on  the  body  shape  and  flight  conditions. 

For  blunt  bodies  at  hypersonic  speeds,  boundary  layer  mixing  and  the 
disturbed  pressure  field  are  the  dominant  influences  in  the  base  region  near 
the  body  and  immediately  downstream,  say  for  distances  up  to  the  order  of 
10  body  diameters.  In  the  far  field  beyond  this,  the  elevated  pressures 
induced  by  the  shock  system  of  the  body  decay  to  the  level  of  the  ambient 
pressure,  and  the  shock-induced  disturbances  which  persist  far  downstream 
are  the  primary  factors.  The  order  of  the  distance  in  which  the  pressure 
downstream  of  a  blunt  axisymmetric  nose  reaches  the  ambient  level  may  be 
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estimated  from  blast-theory,  which  yields  roughly  x/D  ~0.  07M  4 / (p/p,^- 

2 

0.4),  or  for  return  to  ambient  pressure,  x/D  ^  M  / 1 0 .  Therefore,  for 

00 

Mach  numbers  on  the  order  of  20,  downstream  positions  beyond  40  or  50 
body  diameters,  measured  from  the  nose,  may  be  considered  to  constitute 
the  far  field  of  the  wake  and  may  be  treated  as  a  uniform-pressure  field. 
Changes  in  the  remaining  disturbed  properties  can  be  treated  only  by 
considering  diffusive  processes  such  as  those  of  viscosity,  thermal  con¬ 
ductivity  and  mass  diffusion.  Moreover,  under  hypersonic  conditions,  it 
may  be  necessary  to  take  into  account  the  effects  of  chemical  rate  processes. 

Basic  discussions  of  constant-pressure  viscous  flows  of  wake  or  jet 
type  are  given,  for  example,  by  Schlicting  (Reference  2)  for  constant-density, 
and  Pai  (References  3  and  4)  for  compressible  flow.  Some  discussions  or 
analyses  with  more  direct  bearing  on  hypersonic  axisymmetric  wakes  are 
contained  in  the  recent  works  of  Feldman  (Reference  5),  Lew  and  Langelo 
(Reference  6),  Ting  and  Libby  (Reference  7),  M.  and  R.  Goulard  (Reference 
8)  and  Bloom  and  Steiger  (Reference  9).  Laminar  flows  which  are  chemically 
in  equilibrium  or  frozen  are  considered  in  References  5-8.  Ting  and  Libby 
(Reference  7)  also  treat  turbulent  flows.  Bloom  and  Steiger  (Reference  9) 
deal  with  two-dimensional  symmetric  and  asymmetric  cases  as  well  as 
axially  symmetric  flow,  considering  both  laminar  and  turbulent  conditions, 
and  considering  chemical  rate  processes. 
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The  present  work  deals  with  the  behavior  of  the  uniform-pressure 
far  field  wake.  The  flow  is  treated  as  a  continuum,  by  the  boundary-layer 
approximations,  utilizing  integral  method  techniques  according  to  the 
approach  developed|in  Reference  9.  Chemical  rate  processes  of  air  relative 
to  atom-molecule  dissociation  and  to  ionization  are  taken  into  account. 

A  pertinent  general  discussion  of  the  nonequilibrium  thermochemistry 
of  air  has  been  given  by  Bloom  (Reference  10).  Experimental  rate  constants 
for  the  important  processes  cited  by  Bloom  were  those  summarized  by 
Wray,  Teare,  Kivel  and  Hammerling  (Reference  11)  and  by  Teare  and 
Dreiss  (Reference  12). 

Although  the  analysis  to  be  presented  holds  in  general  for  flows  which 
are  non-isoenergetic,  that  is,  of  non-uniform  stagnation  enthalpy,  the  far 
field  wake  may  be  treated  as  approximately  isoenergetic  at  its  initial  station 
since  the  fraction  of  energy  absorbed  upstream  by  the  body  is  negligible. 

In  effect  this  approximation  implies  the  neglect  of  upstream  heat  transfer 
to  the  body  and  of  deviations  of  the  Prandtl  number  from  unity  upstream. 

In  contrast,  wakes  stemming  primarily  from  non-adiabatic  boundary  layers 
over  the  body  or  with  Prandtl  numbers  different  from  unity  must  be  considered 
as  non-isoenergetic  initially. 

Other  approximations  concerning  the  transport  properties  of  the  flow 
will  be  specified  in  the  course  of  the  analysis.  In  particular,  the  approxima¬ 
tions  made  in  dealing  with  mass  diffusion  in  a  multicomponent  mixture  are 
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given  in  Appendix  A. 

The  governing  equations  in  the  wake  are  of  boundary  layer  type  and 
therefore  are  parabolic  in  their  mathematical  structure.  As  a  result,  the 
flow  field  downstream  of  a  given  initial  section,  designated  say  as  x  =  0,  is 
fully  determined  once  the  flow  properties  at  x  *  0  are  specified.  Thus  the 
downstream  flow  depends  upon  the  upstream  flow  only  insofar  as  the  up¬ 
stream  flow  establishes  the  initial  conditions  at  x  ■  0.  It  is  believed  that 
the  wake  properties  will  not  be  strongly  dependent  upon  the  detailed  abera- 
tions  of  the  profiles  of  the  variables  specified  at  x  =  0.  Therefore  the  initial 
profiles  will  be  represented  in  a  smooth  simplified  form  expressible  in  terms 
of  a  polynomial  consistent  with  the  polynomial  profiles  utilized  in  the  inte¬ 
gral  method.  Conditions  furnished  by  estimates  of  the  upstream  flow  are 
matched  at  the  central  axis  and  at  the  outer  edge  of  the  wake.  This  approach 
is  considered  satisfactory  also  because  the  upstream  conditions  are  not 
known  with  accuracy  sufficient  to  warrant  a  more  complex  calculation  at  this 
time.  It  may  be  deemed  desirable  to  make  a  more  accurate  wake  calcula¬ 
tion  when  more  accurate  upstream  properties  become  available. 

The  general  features  of  hypersonic  flow  around  blunt  bodies  with 
regard  to  the  properties  of  the  shock  layer,  shock-shape,  and  near-field 
wake  are  reasonably  well-known,  at  least  for  the  purpose  of  providing  initial 
data  for  the  far-field  wake.  One  approximate  method  for  obtaining  the  re¬ 
quired  initial  data  consists  of  estimating  a  reasonable  shock  shape  for  a 
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given  body,  and  determining  the  conditions  in  the  region  directly  behind  the 
shock,  assuming  that  chemical  equilibrium  is  closely  achieved  near  the 
shock.  The  subsequent  properties  within  stream  tubes  are  then  estimated 
by  neglecting  transport  effects  and  assuming  isentropic  expansion  to  a  given 
pressure,  say  the  ambient  pressure  of  the  far-field  wake.  Isentropic  ex¬ 
pansion  occurs  in  the  limiting  cases  for  which  the  thermodynamic  processes 
are  not  rate-dependent,  that  is,  for  thermodynamic  and  chemical  equilibrium 
or  frozen  conditions.  This  point  of  view  is  supported  by  the  work  of  Bloom 
and  Steiger  (Reference  13)  who  demonstrated  that  the  flow  in  the  inviscid 
supersonic  region  of  hypersonic  blunt  bodies  will  be  frozen  with  respect  to 
atom-molecule  compositipn  for  flight  velocities  between  15,000  and  25,000 
fps  at  altitudes  between  150,  000  and  250,  000  feet. 

Typical  profiles  of  such  initial  data  obtained  by  Ting  and  Libby 
(Reference  7)  for  a  body-size  and  flight  condition  different  from  that  to  be 
considered  here  are  reproduced  in  Figure  2  for  illustration.  Velocity  pro¬ 
files  obtained  for  equilibrium  and  frozen  expansions  are  shown  together  with 
third-degree  polynomial  representations  of  the  initial  data.  In  these  curve- 
fits  the  method  of  least  squares  was  used  to  establish  the  value  of  the  velocity 
at  the  axis  and  to  give  a  value  for  the  radial  extent  of  the  profile,  that  is,  the 
initial  boundary  layer  thickness.  The  initial  profiles  of  atom  and  electron 
concentrations  obtained  by  Ting  and  Libby,  assuming  that  the  concentrations 
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of  species  per  unit  mass  remain  constant  (frozen)  in  the  expansion,  also 
are  reproduced  in  Figure  2  together  with  curve-fits  of  these  profiles. 

Reasoning  is  presented  whereby  the  effect  of  turbulence  in  the  wake 

may  be  estimated.  Of  course,  this  procedure  is  semi-empirical.  It  hinges 

on  the  introduction  of  the  turbulent  eddy  transport  parameters  for  shear,  heat 

flux  and  mass  -  diffusion.  Furthermore,  the  assumption  is  made  that  the 

turbulent  Prandtl  number  and  Lewis  number  are  roughly  equal  in  magnitude 

to  their  counterparts  in  laminar  flow.  Thus,  only  the  eddy  viscosity  \± 

must  be  estimated  explicitly.  In  the  integral  method,  it  develops  that  \x  is 

required  to  be  evaluated  only  at  the  symmetry  axis,  where  it  is  designated 

as  n  .It  has  been  well  established  that  the  turbulent  eddy  viscosity  u.  is 
rto  t 

much  larger  than  its  laminar  counterpart  \x  ,  at  least  in  incompressible  flow. 
In  fact,  the  ratio  ^  may  be  as  large  as  several  hundred  in  incompres¬ 
sible  flow.  There  is  no  reason  to  expect  this  order  of  magnitude  to  change 
radically  in  compressible  flow. 

Based  on  reasoning  of  Mager  (Reference  16)  for  surface  boundary 
layers  and  Ting  and  Libby  (Reference  7)  for  axisymmetric  wakes,  a  somewhat 
rational  estimate  of  the  effect  of  compressibility  on  can  be  made.  In  the 
wake  case,  experimental  checks  of  this  procedure  have  yet  to  be  made;  there¬ 
fore  the  procedure  must  be  considered  arbitrary  at  this  time.  However,  it  is 
reasonable  and  thus  warrants  consideration.  In  any  event,  this  estimate  of  the 
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effect  of  compressibility  on  p^  does  not  show  a  change  in  the  order  of  magni¬ 
tude  of  /  p  from  that  of  incompressible  flow. 

It  will  be  seen  that  the  effect  of  turbulence  appears  in  two  major 
ways.  First,  the  streamwise  length  of  the  wake  is  inversely  proportional  to 
the  molecular  viscosity  p  at  the  axis  in  laminar  flow,  and  is  inversely 
proportional  to  the  eddy  viscosity  ptQ  at  the  axis  in  turbulent  flow.  Therefore 
the  streamwise  extent  of  the  wake  with  respect  to  the  velocity  field,  state 
variables  and  concentration  of  major  components,  is  diminished  roughly  in 
the  ratio  p  /  p^  due  to  turbulence.  That  is,  the  turbulent  wake  may  be 
several  hundred  times  shorter  than  the  comparable  laminar  wake. 

Secondly,  in  the  equation  governing  species  concentrations,  the 

ratio  of  the  chemical  production  rate  term  to  the  diffusion  term  is  inversely 

proportional  to  p  in  laminar  flow  and  inversely  proportional  to  p  in  turbulent 
o  to 

flow.  Therefore,  the  importance  of  chemical  rate  processes  relative  to  mass 
diffusion  is  diminished  by  orders  of  magnitude  in  turbulent  flow. 

For  the  altitudes  (150,000  -  250,000  ft.)  and  flight  velocities 
(18,900  -  20,400  fps)  of  direct  interest  here,  it  will  be  seen  that  the  chemical 
rate  effects  are  negligible  with  respect  to  mas s -diffusion  effects  even  in 
laminar  flow  for  the  primary  air  components,  but  not  for  the  electron-producing 
process.  That  is,  the  flow  will  be  chemically  frozen  for  the  main  air  species 
under  these  conditions,  although  it  will  not  be  frozen  in  this  manner  at  suffi¬ 


ciently  low  altitudes. 
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The  necessity  of  using  a  chemical  rate  process  for  calculating 
electron  densities  in  a  hypersonic  flow  field  is  fully  brought  out  by  the  pre¬ 
sent  investigation.  The  results  indicate  that  the  ion  species  (  a^jQ+  ) 
definitely  in  a  nonequilibrium  state  for  both  the  supersonic  expansion  region 
and  the  wake,  for  flight  velocities  of  about  20,000  fps  and  altitudes  between 
150,  000  and  250,  000  feet. 

Very  recently  Rosenbaum  and  Bloom  (Reference  17)  completed  an 
investigation  of  the  inviscid  nonequilibrium  molecular  dissociation  and 
ionization  flow  about  a  body  with  a  6  inch  hemispherical  nose  radius  for  free 
stream  flight  conditions  of  200,000  feet  and  20,200  fps.  Their  results  show 
that  for  these  flight  conditions  (1)  all  species  (  a^,  a^,  a^Q+)  attain 
equilibrium  in  a  relatively  close  distance  behind  the  bow  shock,  (2)  the  main 
species  (a^>  a^)  are  frozen  in  the  supersonic  expansion  region  and  (3)  the 
ion  species  ( cx a  nonequilibrium  state  which  remains  closer  to 
equilibrium  than  the  frozen  values .  The  eLectron  distribution  about  the  body 
is  presented  in  Figure  16. 
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II  ANALYSIS 

The  following  system  of  equations  is  assumed  to  govern  the  flow 
in  an  axiaLLy- symmetric  Laminar  wake,  shown  schematicaLLy  in  Figure  1. 
The  assumptions  utiLized  in  representing  mass  diffusion  in  the  muLticom- 
ponent  reacting  air  mixture,  are  discussed  in  more  detaiL  in  Appendix  A. 


OveraLL  mass  conservation: 


(puy)  +  ( pvy)  =  0 

x  y 

(1) 

Momentum: 

(pu  2  y)  +  (puvy)  = 

x  y 

(PyVy 

(2a) 

p  =0 

y 

(2b) 

Energy: 

(puyH)x  +  ( pvyH) y  =  £(p/P)yHyJ 

|y  +  |(P-l)(|i/P)yuuy 

+ 

y 

+  [(Le-l)(p/P)y  hiaiyJy  (3) 

i 

Mass  conservation  of  species  i: 


(puya^  +  (pvya^y 


(Lep/P)yaiy  J  y  +  pywA 


(4) 


Notation  not  defined  in  the  text  is  given  in  an  attached  List  . 

Before  the  supplementary  thermodynamic  and  chemical  relations 
are  discussed  the  boundary  conditions  wilL  be  stated  and  some  resuLts  not 

dependent  upon  thermochemistry  derived.  The  minimum  number  of 
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independent  boundary  conditions  is  shown  enclosed  in  square  brackets. 
The  remaining  boundary  conditions  are  derived  from  the  differential 
equations  or  are  prescribed  on  physical  grounds  in  the  usual  manner  of 
the  integral  method. 

Boundary  Conditions: 

at  y  =  0:  [v  =  o]  .  £y  =  o]  ,  jHy  =  o]  ,  |«iy  =  §  (5a) 


(pii/2)  ux  =  (jiuy) 


(5  b) 


at  y-~>6: 


(pu/2)Hx  =  [WP)Hy]y  +  QP-l)(pt/P)uuyJ 

+  |~(Le-l)(,i/P)  ^hiaiyly  (5c) 

i  -J 

(pu/i)oix  =  flLep/Plaj,!  ^  +  pw./2  (5d) 

u  s  uq(x)  ,  H  i  Ho(x)  ,  ai  i  a  io(x)  (5e) 

[“  =  ue].[H  *  He].  [aj  =  aie]  (5f) 

“y  "  uyy  =  0  ’  Hv  =  V  =  0  '  “iy  =  aiyy  =  0 


Integral  Equations: 

The  integral  equations  following  from  (1)  -  (5)  are: 

r&  -- 

Qx  =  0  or  0  =  0  (constant)  ;  0  =  J Q  pu(l-u)y  dy  (6a) 

(0E)X  =  0  or  eE  =  ©Ec( constant);  ©E  =  J ^  pu(  1-H)y  dy  (6b) 
(0i)x  =  U/ue)  f  pw.y  dy  ;  0.  =  f  pu  (a.g-a.)  y  dy  (6c) 
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The  foLLowing  transformations  are  now  introduced  into 


the  integral,  equations  (6)  and  boundary  conditions  (5): 

Z 

O 


mdm  =  py  dy  ,  m  =  &mN  ;  6  /  2  =  J  m(  1  /  p)  mdm  ( 7) 

ds  =  ;ro  dx/(peueecAie)  ,  s-sc  =  (  f  ho  dx)/(peueoc/jie) 


It  is  observed  that  0,  0  and  0.  have  the  dimension  of 

E  i 

(Length)2  ;  m,  6,  and  have  the  dimension  of  (Length);  whiLe  s  and  N 


are  dimensionLes  s . 


The  foLLowing  working  forms  are  derived  from  (6)  and 


(  5b-d): 


6  2  f1  u(l-u)  NdN  =  9, 
m  Jo  < 


(6  /  0  )  u  u  =  2vL11(.. 

'  rrv  c  o  os  NNo 


(  8a) 

(  8b) 
(9a) 


6  f  u(l  -H)  NdN  =  0 
n \Jo  EC 

<6j/Sc)GoHoS  M2/Po)SNNo  +  (4/Po,(Po-1)(ue  /2He)%“NNo 

+  l2/PJU-.0-l)^  h.c  aiNNo 


0.  =  (p  u  0„/u  )(l/u.  )(6  ‘/u  )  f1  pw.NdN 

is  'Ke  e  c m'  e'  Jo  1 


(9b) 


(10a) 


(6  /  0  )  u  a . 

'  m  '  c  o  ios 


(  2Leo/P0)  aiNNo  +  (  l/p  0  ^MPe^c-V.)  • 


O  e 


(6  2  /0  )  w. 

'  m  '  c'  k 


(10b) 
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Direct  consideration  of  the  original  differential  equations  shows 


that  for  the  approximation  Le  =  1 ,  the  energy  equation  (3)  becomes 
independent  of  concentrations  except  indirectly  through  jjl  and  P.  This 
simplification  carries  over  to  the  energy  boundary  condition  (9b),  where, 
in  addition,  the  expLicit  appearance  of  (ji  has  been  eliminated  by  the 
transformation  of  x.  For  the  further  approximation  P  =  1,  the  Crocco 
integral  relation  H  =  A  +  Bu,  where  A  and  B  are  constants,  is  applicable 
in  view  of  the  axial  symmetry.  By  satisfying  the  relations  H  =  H  and 


oc 


H  when  u  =  u  and  u  ,  this  relation  may  be  expressed  in  the  form: 


e 


oc 


e 


H  -  H  u 


e 


u 


e 


(11) 


Likewise,  when  the  effects  of  chemical  rates  are  negligible  and 


the  term  involving  w.  in  (4)  is  dropped,  and  moreover  Le/  P  =  1,  the 
following  integraL  of  (4)  is  derived: 


u  -  u 
e 


(12) 


This  relation  is  comparabLe  to  (11)  and  satisfies  the  conditions 


and  a  when  u  =  u  and  u  . 
ie  oc  e 
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It  is  recaLLed  that  the  integral,  method  for  soLving  partial  differen¬ 
tial  equations  is  an  approximate  method  which  is  used  in  cases  where  more 
exact  calculations  are  difficult  to  obtain,  or  where  a  moderate  sacrifice 
in  accuracy  is  permissible  in  the  interest  of  simplicity.  It  involves  the 
assumption  of  functional  form  s  of  the  dependent  variables.  These  forms 
are  chosen  in  a  reasonable  fashion  to  satisfy  a  number  of  boundary  condi¬ 
tions,  but  are  permitted  to  retain  undetermined  parameters,  which  are  then 
utilized  to  satisfy  the  integral  equations  and/or  additional  boundary  condi¬ 
tions,  or  other  conditions  derivable  from  the  differential  equations.  These 
other  conditions  may  be  derived,  for  example,  by  requiring  that  the 
differential  equations  be  satisfied  at  certain  interior  points  or  regions 
away  from  the  boundaries.  On  the  other  hand  they  may  take  the  form  of 
additional  integral  equations  derived  by  piecewise  integration  of  the  differen¬ 
tial  equations  over  one  of  the  variables  (in  boundary- laye r  problems  usually 
the  coordinate  normal  to  the  main  stream),  or  by  applying  weighing  factors 
to  the  original  differential  equations  prior  to  integration.  The  weighing 
functions  may  be  powers  of  the  normal  coordinate  or  of  the  streamwise 
velocity,  for  example.  Clearly,  the  larger  the  number  of  independent 
conditions  which  are  satisfied  by  the  assumed  functions,  the  greater  will 
be  the  accuracy  of  the  solution.  A  compromise  between  simplicity  and 
accuracy  influences  the  choice  of  the  specific  procedure  to  be  employed  in 
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a  given  problem;  this  choice  usuaLLy  invoLves  a  degree  of  judgment  based 
upon  experience.  More  accurate  test  caLcuLations,  limiting  cases,  special 
cases  whose  mathematically  exact  solutions  are  known,  and  experimental 
results  provide  yardsticks  for  evaluating  the  accuracy  or  reasonableness 
of  a  solution  obtained  by  the  integral  method. 

In  the  present  analysis  it  is  assumed  that  the  overall  behavior  of 
the  profiles  of  velocity,  stagnation  enthalpy  and  species  concentration  can 
be  depicted  adequately  by  fourth-degree  polynomials  in  the  normalized, 
transformed  radial  coordinate  N,  and  that  only  one  undetermined  parameter 
included  explicitly  in  each  profile  is  sufficient  to  represent  the  effect  of 
streamwise  variations.  Another  undetermined  parameter,  the  wake- 
thickness  6,  assumed  to  be  the  same  for  all  flow  variables,  is  included 
implicitly  through  the  definition  of  N. 

This  selection  of  the  profile-forms  implies  that  the  initiaL  data 
for  each  variable  can  be  represented  with  sufficient  accuracy  by  a  fourth- 
degree  polynomial  with  one  constant,  aside  from  the  wake -thickne s s 
common  to  all  variables,  available  for  adjustment.  In  practice,  for  a 
symmetric  wake,  the  adjustable  quantity  in  each  profile  will  be  the  value 
of  the  variable  evaluated  at  the  center  axis.  For  additional  accuracy  in 
fitting  the  initial  data  or  in  representing  the  flow  field,  additional  para¬ 
meters  may  be  included  in  the  profile  representation.  Beyond  the  initial 
station  these  may  be  evaluated  in  the  manner  outlined  above. 
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Assumed  Profil.es: 

The  assumed  profiles,  which  satisfy  boundary  conditions  (5a),  (5f) 
and  definitions  (5e)  are  seen  to  be  the  same  for  u,  H  and  a^>  namely: 


u  -  u  H  -  H 
o  o 


a  .  -  a  . 

_ 1  IQ 

u-u  H-H  a-  -  a ;  ^ 
e  o  e  o  le  io 


or  in  another  form, 


u  -u  H  -  H  a  -  a. 
e  e  ie  i 


2  3  4 

6 N  -  8N  +  3N 


ue  ~  uo  ^e  ~  H0  aie“  aio 


2  3  4 
1  -  6N  +  8N  -  3N 


(13a) 


(13b) 


The  profile  forms  (13a)  are  plotted  in  Figure  3. 

The  relations  in  (13)  should  not  be  confused  with  the  special  inte¬ 
grals  (11)  and  (12).  However,  they  do  lead  to  the  following  relations  based 


on  the  definitions  of  0,  0  and  0.  given  in  (6): 

E  i 


0  = 


1  ~  go 

1  -  H 


E 


l  -  a. 


k  a  .  -  a  .  . 

ie  ioj 


0. 


(14) 


Now  if  we  choose  to  satisfy  the  momentum  and  energy  integral 


equations  (8a)  and  (9a),  which  state  that  both  0  and  0  are  constants  along 

E 

the  flow,  we  derive: 


1  -  u 

o 

1  -  H 
o 


0 

<T 


u 


oc 


'Ec 


1  -  H. 


oc 


(15) 
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which,  with  assumption  (13b),  yieLds  (11).  Thus  it  is  seen  that  the  assum¬ 
ption  that  the  veLocity  and  stagnation  enthaLpy  profii.es  are  of  the  same  form, 
and  the  satisfaction  of  the  momentum  and  energy  integral  relations,  lead 
to  the  special  relation  (11)  ,  but  without  any  assumptions  concerning  the 
transport  properties  Le  and  P. 

Likewise,  for  the  case  of  negligible  chemical  rate  effects,  the 
assumption  that  the  velocity  and  species  concentration  profiles  are  of  the 
same  form,  and  the  satisfaction  of  the  momentum  and  concentration 
integral  relations  yields  (12)  but  without  any  assumptions  concerning  Le/P. 

Viewed  in  another  way,  it  is  seen  that  the  approximate  relations 
assumed  between  velocity,  enthalpy  and  concentrations  become  exact  in 
the  special  cases  cited. 

The  polynomial  in  N  assumed  in  (13a)  varies  monotonicaLly  from 
zero  at  the  axis  to  unity  at  the  edge  of  the  wake.  This  is  consistent  with 
the  initial  conditions  with  which  we  wish  to  deal  here.  As  a  consequence 
of  this  assumption  it  is  seen  from  (6)  that  if  uq  =  u^  or  Hq  =  at  any  one 
streamwise  station,  it  is  required  that  0  =  0  or  0  =  0  throughout  the  flow. 

hj 

The  condition  0  =  0  is  satisfied  onLy  in  the  triviaL  case  of  uniform  flow. 
However,  0g  =  0  implies  a  flow  with  uniform  stagnation  enthalpy,  which  is 
quite  possible  in  view  of  the  adiabatic  condition  of  the  wake  with  respect 
to  its  surroundings.  Analogous  considerations  hold  for  the  distributions 
of  species  -  concentration  when  chemical  rate  effects  are  negligible. 
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Integral,  method  solutions:  veLocity  and  wake-thickness  parameters 
We  choose  to  satisfy  the  momentum  integral,  equation  (8a)  and  the 
momentum  condition  aLong  the  axis  (8b)  in  order  to  determine  the  trans¬ 
formed  wake  thickness  6^  and  the  veLocity  parameter  uq  (s).  Substitution 
of  the  assumed  profiLe  (13)  into  (8)  yields: 


0/6  2  =  (1  -  u  )(10  +  11  u  )/210 

c'  m  o  o 

(16) 

U  u  =  (  4/35)(l  -  u  )2  (10  +  11  a  ) 
o  os  o  o 

(17a) 

A  pLot  of  (16)  is  shown  in  Figure  4. 

Equation  (17a)  is  separable  and  its  integration  yieLds  an  equation 

of  the  form: 

s  =  f(uQ)  +  constant  (17b) 

where 

f(uo)  r  (5/12H 1/1  -  uo)  -  (25/126)  in  £(10  +  1 1  uq)(  1  -  uq)J  (17c) 

Since  the  absolute  origin  of  the  axiaL  coordinate  x  or  s  is  arbitrary 

within  an  additive  constant,  we  may  define  the  position  x  =0  or  s  =0  to  be 

the  Location  at  which  a  given  initial  value  of  uq  exists.  With  this  initiaL 

value  denoted  as  uQC>  (1Tb)  can  be  expressed  as 

s  =  f  (  u  )  -  f  (  u  )  (1  7d) 

o  oc 

Equation  (L7d)  is  pLotted  in  Figures  for  uqc  =■  0.  This  curve  is 
generaL  in  the  sense  that  if  we  consider  a  particular  value,  say  uQa,  which 
corresponds  to  a  coordinate  vaLue  s  ,  the  curve  gives  the  subsequent  or 
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prior  values  of  u  which  correspond  to  coordinate  increments  s  -  s  , 
o  r  a 

regardless  of  the  absolute  location  of  s^.  As  s-  >oo,  u— ~->l. 

Stagnation  enthalpy  parameter: 

Now  if  we  choose  to  determine  H  (s)  by  satisfying  the  energy 

integral  equation  (9a)  we  derive  (15)  as  previously  mentioned.  That  is 

(1-  H  )/(l-  H  )  =  (l-u  )/(l-u  )  (18) 

o  oc  o  oc 

where  H  may  be  specified  independently.  Of  course,  when  Hqc  =  1  it 
is  required  that  Hq  =  1  in  accordance  with  the  previous  discussion  of  flow 
with  uniform  H.  This  solution  for  Hq,  derived  from  the  integral  equation 
which  essentially  expresses  an  average  across  the  flow,  does  not  depend 
upon  P  or  Le .  On  the  other  hand,  an  idea  of  the  Prandtl  number  -  Lewis 
number  influence  can  be  obtained  by  alternately  satisfying  the  energy 
differential  equation  along  the  center  axis.  For  this  purpose  (9b)  may  be 
combined  with  (16),  (17a)  and  the  results  of  (13)  to  yield  the  following 
expression: 


dH 

o 

i 

^No 

2(P-  1) 

o 

2 

U 

e  

(Le  -  1) 
o 

Ik  . 

i  io  iNNo 

du 

o 

Po 

+ 

12(1-  a  ) 

o' 

P 

o 

uo  T 

2H  P 

e  o 

o 

1 

1 

U-H0> 

2(P-1) 

2 

u 

e 

(Le  -1) 
o 

h.  ( a  .  a  .  ) 

i  io  le  io 

1  &° 

1 

<*-%>  + 

P 

o 

uo 

2H  P 

e  o 

d-do) 

(19a) 


(19b) 


1-  H 


1-  u. 


1-P 


1  -  H  u 
o  e 


u 


l-un  H„ 
o  e 


z 


(Le  -1)  4—  h  (  a.  -  a  .  )  (19c) 
o  1  io  ie  10 


1-  u 
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For  P  =1  and  Le  =1,  this  equation  yields  (18).  A  rough  estimate 
o  o 

of  the  effects  of  deviations  of  and  from  unity  can  be  obtained  by 

making  an  a  priori  estimate  of  the  order  of  magnitude  of  the  two  additional 
terms  on  the  right  hand  side  of  (19c),  and  then  integrating  (19c)  as  a  linear 
equation.  This  is  equivalent  to  forming  an  integral  equation  from  (19c)  and 
obtaining  a  first  iterative  solution.  The  procedure  is  shown  in  Appendix  A. 

The  considerations  in  Appendix  A  show  that,  for  constant  average 
values  of  the  second  and  third  right  hand  terms  in  (19c),  whose  sum  is 
denoted  by  B,  (19c)  may  be  integrated  to  yield 

1  -  u 


1  -  H  1  -  H 
_ 9.  r  _  QS- 


1  -  u 


B  in 


oc 


1  -  u 


oc 


1  -  u 


(20) 


where  B  is  a  negative  quantity  whose  magnitude  is  on  the  order  of 

(1  -  P  )  /  PQ  for  hypersonic  conditions  in  the  far-field  wake  and 

Le  -  1  ~  1  -  p  . 
o  o 

Representative  values  of  H  as  a  function  of  u  ,  calculated  from  (20), 
are  shown  in  Figure  6  for  B  =0,  corresponding  to  P  =  Le  =1,  and  the 


typical  value  J B  J  =  0.3.  Cases  are  shown  for  an  initial  stagnation  enthalpy 

equal  to  that  of  the  free  stream,  i.  e.  ,  H  =1,  and  for  a  small  initial 

oc 

deviation  in  H  ,  i.  e.  H  =1-0.  2(1  -  u  ).  The  effect  of  P  /  1  and 


oc 


oc 


Le  /  1  is  to  encourage  decreases  in  H  below  H  . 
o  o  e 

additional  decreases  is  about  5  per  cent  of  H  . 

e 


The  order  of  these 
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Species  concentration  for  negligible  chemical  rates  (frozen)  : 

It  is  possibLe  to  show  that  for  a  cLass  of  technically-important 
wake  flows  including  those  of  direct  interest  here,  certain  chemical 
processes  occur  with  sufficient  sLowness  that  the  rate  effects  are  negligible 
with  respect  to  the  corresponding  diffusion  effects.  Such  flows  are  termed 
"frozen*1  with  respect  to  these  particular  processes.  Specifically,  the 
processes  involving  the  dissociation  of  02  and  N2  and  recombination  into 
02  and  N2  will  be  frozen  at  the  altitudes  and  flow  velocities  of  speciaL 
interest  here.  This  is  not  necessarily  the  case  for  trace  components  such 
as  NO  or  for  the  electron-producing  processes.  However,  the  02  and  N2 
processes  are  particularly  important  since  they  dominate  the  state  variables, 
such  as  density  and  temperature,  of  the  air  mixture. 

Solutions  for  the  concentrations  of  the  frozen  species  under  the 
influence  of  mass  diffusion  can  be  pursued  without  a  discussion  of  the 
thermochemical  processes  of  air  components.  Therefore  the  discussion 
of  the  thermochemistry  wiLL  be  given  Later. 

With  the  w.  terms  negLected,  (10a)  and  (  10b)  take  the  form: 

0.  =  0.  (constant)  (21a) 

1  1C 

2 

(6  /  0  )  u  Q  =  (  2  Le  /  P  )  a 

m  c  o  ios  o  o  iNNo 


(21b) 
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The  solution  of  these  equations  bears  a  strong  resemblance  to 
those  governing  the  stagnation  enthalpy  parameter  Hq  and  discussed 
previously.  For  example,  when  Le  /  P  =1,  (12)  follows  from  the  differen¬ 
tial  equations.  Likewise,  with  the  assumed  profiles  (13),  satisfaction  of 
the  integral  equation  (21a)  above  yields  (12)  with  no  specification  of  Le/P. 
Furthermore,  the  alternate  satisfaction  of  the  differential  equation  along  the 
axis  may  be  accomplished  by  the  use  of  (21b)  coupled  with  (16),  (17a)  and 
(13),  which  yields 

da  .  Le  ( a .  -  a .  ) 

io  o  ie  lo  . 

-  - - -  (21c) 

dSo  Po  U’V 

Integration  of  (  21  c),  assuming  Le  /P  =  constant,  yields 


and  clearly  shows  the  order  of  the  effect  of  Le  /P  ^  1.  Since  Le  /p  ^2 

o  o  o'  o 

in  high  temperature  air,  its  effect  on  the  diffusion -c  ontr oiled  species 

concentrations  can  be  appreciable.  Eq.  (22)  is  shown  plotted  in  Figure 

7a  for  Le  /P  =1  and  2.  Corresponding  pLots  in  terms  of  the  streamwise 
o  o 

coordinate  are  shown  in  Figure  7b. 
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Species  concentration  for  very  rapid  chemical  rates  (equilibrium): 

In  the  Limiting  case  in  which  the  rate  parameters  of  the  chemical 
processes  are  very  Large,  the  species  concentrations  are  uniquely  deter¬ 
mined  by  the  usual  relations  governing  equilibrium  thermodynamics 
(mass-action  Law  for  each  reaction;  equation  of  state;  relation  between 
enthalpy,  concentrations,  and  temperature;  conservation  of  atoms  in 
the  course  of  the  reaction)  once  the  pressure  and  enthalpy  are  specified. 

In  this  case  the  quantity  w.,  which  represents  the  local  rate  of  change 
of  concentration  of  a  given  species  due  to  chemical  change  alone  ,  becomes 
an' indeterminate  form  formed  by  the  product  of  a  very  large  rate  parameter 
and  the  difference  between  the  forward  and  reverse  proce  ss  terms  which 

approaches  zero.  There  is  no  need  to  evaluate  w  directly,  since  in  this 

i 

case,  the  species  conservation  equation  (4)  simply  expresses  the  fact 

that  if  we  consider  a  mass  element  in  transit,  the  change  in  species 

concentrations  within  the  element  is  due  to  the  chemical  change  effect, 

signified  by  w  ,  plus  the  diffusion  effect.  Therefore,  when  the  concen- 
i 

trations  are  known,  as  in  the  case  of  chemical  equilibrium,  the  convective 
term  and  diffusion  term  in  (4)  may  be  evaluated,  and  thus  w.  may  be 
determined  from  (4)  if  necessary.  In  practice  (4)  may  be  disregarded 
in  the  case  of  chemical  equilibrium,  since,  in  effect  it  is  automatically 


satisfied. 
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Thermochemistry  of  air: 

At  the  energy  Levels  involved  in  the  wake  flows  under  consideration 
here,  the  state  variables  such  as  density  and  temperature  of  the  air  mixture, 
are  dominated  by  the  primary  components  Oz  ,  N2  O  and  N.  ALL  other 
components  are  expected  to  be  present  only  as  traces  whose  energy  changes 
are  negligibly  small  with  respect  to  those  cited.  Of  course  the  trace  com¬ 
ponents  are  important  from  the  standpoint  of  their  influence  on  radiation  or 
electron  densities.  However,  the  primary  processes  involving  Q2~-  v  20 
and  N2  v  —  2N  may  be  considered  first,  by  ignoring  the  trace  components. 
Then  the  trace  processes  can  be  investigated  utilizing  the  previously-deter¬ 
mined  state  and  concentrations  of  major  species  as  input  information. 

In  the  relations  to  be  given,  the  individual  species  of  the  air  mixture 
are  assumed  to  behave  as  perfect  gases.  Furthermore  the  translational 
rotational  and  vibrational  degrees  of  freedom  of  the  moLecules  are  assumed 
to  be  in  equilibrium,  despite  the  fact  that  dissociation  and  ionization  pro¬ 
cesses  may  be  rate -dependent . 

The  principles  of  chemical  kinetics  and  statistical  thermodynamics 

are  applied  to  the  following  reactions  of  three -body  collision  type: 

kdN 

N2  +  X~_Z±2N  +  X  (23a) 

k 

rN 
k 

dO 

02  +  X'~~ — ^  20  +  X 

k 

rO 


(23b) 
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where  X  signifies  the  cataLyst,  which  here  is  assumed  to  consist  of  any 
other  particLe  contained  in  the  mixture.  Refinements  of  the  type  which 
consider  the  influence  of  various  types  of  cataLysts  on  the  reaction  rates 
are  not  considered  here.  This  wouLd  require  consideration  of  a  separate 
reaction  and  rate  process  for  each  catalyst,  and  wouLd  resuLt  in  coupling 
which  wouLd  complicate  the  computational  procedure  considerably. 

The  following  relations  corresponding  to  (23)  are  derived: 

ChemicaL  production  rates: 


w 


=  (kr0Pd0/M2c)(M0! 


-Dn/RT 


a  ^  e 

o2 


-  (p/pdo)  ao 


(24a) 


w 

N 


=  (kMP  /M) 

rN  dN  N  N2 


-dn/RT  /  /  V 

e  1N  -  (p/p,  ) 

dN 


(24b) 


w  =  -  w 

o2  O 


w_T  =  -  w  (24  c,  d) 

N2  N 


where 


l/M  =  a./M. 


(24e) 


Pdi  =  (m./2)(kd./kri)e 


D./RT 


3/2/_R  -Ti/Tx/Evz  ,  E 


=  (m./hhUm.kT)  (t*/t)(i  -  e  1  )(**)*/*. 


(24f) 


1 


O  or  N 
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where  k  and  k  are  the  forward  and  reverse  reaction  rate  parameters 
di  ri 

shown  in  (23),  M.  denotes  moLecuLar  mass  in  mass/mole,  m.  denotes 
1  i 

particle  mass  in  mass/particle,  D.  is  the  dissociation  energy  per  mole 

associated  with  the  atomic  species,  h^  is  PLanck*s  constant,  k  the  gas 

R 

constant  per  molecule,  R  the  gas  constant  per  mole,  T  the  characteristic 

i 

V 

rotational  temperature  of  a  species,  the  characteristic  vibrational 

E 

temperature,  f  the  electronic  partition  function  of  atomic  species,  and 
i 

E 

f  the  electronic  partition  function  for  corresponding  molecular  species. 

k 

Enthalpy  of  the  mixture: 


h  = 


(RT/2M  J  j  3.88  +  0.232  A  +  0.875AXT  +  a  (1.50  -A  )  + 

'  o  (  O2  H  o  o2 


+  aN(1.14)(1.50  -A  ) 


1 


J 


+  a  D/2M  +  a  D  /  2  M 
OO  O  N  N  N 


(25) 


where 


V  , 

V  T.  /T 

-A.  =(T.  /  T)  /  (e  1  -  1)  ;  i  =  Oz  ,  N2 

Of  course  the  number  density  of  electrons  is  equal  to  the  number 

+ 

e  NO 

density  of  ions  in  the  present  approximation.  Therefore  n  =  n  = 

(pa  +)  /  (m  +),  where  m  is  the  particle  mass  which  is  equal  to  M. 

NO  ' 9  x  NO  1  i  1 

divided  by  Avogadro*s  number  and  n.  has  the  units  particles  per  unit  volume. 
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Equation  of  state: 

p  =  pRT/M  (26) 

Numerical  vaLues  for  use  in  the  above  thermochemical  relations 
are  given  inReference  LO  unLess  specified  eLsewhere. 

Electron  production: 

In  view  of  our  interest  in  the  eLectron  density  distribution  in  the 
fLow  fieLd,  we  shaLl  consider  the  following  e  Lectron-producing  reaction  which 
has  been  found  to  be  of  importance  in  experiments  summarized  in  References 
11  and  12: 

N  +  O  ^ 


kdNO+  v  + 

- ==  NO  +  e 


:NO+ 


(27) 


This  leads  to 


WNO+  ~  *krNO+  *  P 


^MNO+  /  MNM0^kdNO+//krNO+  )  aNao"  aNO+ 


(28a) 


where,  with  T  in  R,  and  the  ionization  energy  per  mole: 

V2 

kdNO+//krNO+  =  (1/2)(mNO+  Me~/M0MN^ 


-320/T  -4LL/T  -589/T 

(2  +2e  )  /  (5  +  3e  +  e  ) 


4100/T  4100/T 

(T/8.83)(e  )  /  (e  -1) 


-  (I  ,  -  D  )  /  RT 
NO+  NO+  ' 


(28b) 
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ALL  the  information  required  for  making  calculations  of  species 
concentrations  of  O,  N,  02  ,  N2  and  subsequently  NO+  and  e"  is  now  at 
hand.  For  the  calculations  to  be  presented,  onLy  the  ordinary  differential 
equation  obtained  from  the  partial  differential  equation  along  the  symmetry 
axis  will  be  employed.  The  integral  equations  for  stagnation  enthalpy  and 
species  concentration  wiLL  be  disregarded.  This  choice  is  made  because  the 
flow  field  near  the  axis  is  believed  to  be  of  major  importance,  because  the 
condition  along  the  axis  retains  a  considerable  amount  of  mathematical 
exactitude  and  because  Prandtl  number  and  Lewis  number  effects  are  re¬ 
tained.  In  effect  then  we  are  simpLy  satisfying  the  differential  equation 
along  the  axis,  estimating  that  the  second  derivatives  of  the  variables  are 
proportional  to  the  difference  between  the  outer  vaLues  and  inner  values  of 
the  variables,  and  then  using  the  assumed  profiles  to  interpolate  between 
the  vaLues  at  the  axis  and  the  outer  vaLues. 

Transport  properties: 

The  means  for  evaluating  the  transport  properties  of  gas  mixtures 
are  now  well  known,  particularly  in  the  case  of  binary  mixtures.  Detailed 
discussions  of  these  calculations  are  given  in  References  14  and  15.  The 
expressions  for  the  transport  properties  wiLL  not  be  repeated  here  for  the 


following  reasons: 
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The  thermal  conductivity  X. ,  as  well  as  the  mean  specific  heat  c  , 

P 

a  thermodynamic  property,  appear  together  with  the  viscosity  coefficient 

p  in  the  Prandtl  number  P=  pc  /  \  whereas  the  mean  diffusion  coefficient 

P 

(D..)  ,  which  is  used  here,  appears  in  the  Lewis  number  Le  =  p(D  )  c  /\  . 

av  ..  ij  av  p 

The  parameters  P  and  Le  are  not  evaluated  in  detail  here,  since  it  is  well 
known  that  they  do  not  exhibit  a  strong  dependence  on  the  state  (including  the 
precise  composition)  of  air  or  of  mixtures  of  its  main  components.  That  is, 
for  such  mixtures,  P  takes  on  values  roughly  between  0.7  and  1.0,  while 
Le  assumes  values  on  the  order  of  1.4.  Along  the  lines  previously  discussed, 
constant  average  values  of  P  and  Le  are  assumed  for  the  present  calculations. 

On  the  other  hand,  the  viscosity  coefficient  appears  explicitly.  It 
has  been  determined  that  for  air  in  equilibrium  at  high  temperatures,  the 
viscosity  p  is  correlated  within  an  accuracy  of  10%  by  the  well-known 
Sutherland  relation,  which  can  be  written  in  the  following  form: 


p  =  T  2 


1  +  C 


(29) 


1  +  C/  T( 

where  p  =p/pe>  T  =  T/Te  and  C  =  200OR/T  .  This  relation  is  used  in 
the  present  calculations. 

Turbulent  flow: 

Methods  for  treating  the  turbulent  fLow  of  wakes  and  jets  have  received 


considerable  treatment  in  the  literature  (for  example,  see  Reference  2, 
Chapter  23).  In  effect,  in  turbulent  flow,  each  of  the  molecular  transport 


GENERAL  APPLIED  SCIENCE  LABORATORIES,  INC. 


page  ^9 


effects  such  as  shear  stress,  heat  flux  and  mass  diffusion  velocity,  is 
augmented  by  an  additional  quantity  due  to  the  turbulent  transports  associated 
with  the  turbulent  fluctuations.  Based  on  physical  reasoning  albeit  with  a 
degree  of  arbitrariness,  definitions  can  be  introduced  to  express  the  turbu¬ 
lent  shear,  turbulent  heat  flux  and  turbulent  mass  diffusion  in  forms  analo¬ 
gous  to  the  molecular  transports.  That  is,  the  shear  can  be  set  proportional 
to  the  gradient  of  the  mean  velocity,  the  heat  flux  proportional  to  the  gradient 
of  mean  temperature  and  the  mas s -diffusion  velocity  proportional  to  the 
gradient  of  the  mean  concentration.  The  proportionality  factors  are  termed 

eddy-viscosity  p  ,  eddy-conductivity  \  and  eddy-diffusivity  (D  )  respec- 
t  t  ij  t 

tively,  and  are  determined  by  combinations  of  physical  reasoning  and 
experimental  observation.  The  turbulent  transports  outweigh  the  molecular 
transports  in  importance  in  flow  regions  which  are  not  extremely  close  to 
solid  surfaces,  where  laminarization  takes  place.  In  fact,  in  free  turbulent 
flows,  such  as  wakes  and  jets,  the  molecular  transports  can  usually  be 
neglected  with  respect  to  the  turbulent  contributions  throughout  the  flow. 

In  practice,  the  definition  of  the  eddy  parameters  simply  means 
that  in  turbulent  flows  we  may  utilize  the  laminar  flow  equations,  with  the 
interpretation  that  the  flow  variables  signify  the  corresponding  mean  quanti¬ 
ties  of  turbulent  flow,  and  that  the  transport  parameters  (p,  D  )  in  the 

ij 

laminar  equations  actually  denote  the  sums  of  the  laminar  quantities  and  the 
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corresponding  empirical  turbulent  eddy-parameters  p.  ,  \  and  (D..)  . 

t  t  ij  t 

That  is,  p.  is  replaced  by  (p.  +  p.^)  etc.  in  the  equations  stated  in  (1)  to  (4). 
This  can  be  done  since  the  differential  equations  for  the  mean  turbulent 
motion  are  otherwise  identical  to  those  for  laminar  flow  provided  that  the 
density  fluctuations  are  neglected.  Asa  refinement  the  density  fluctuations 
can  be  retained  but  this  is  not  done  usually;  their  significance  has  not  been 
established . 

An  idea  of  the  magnitude  of  the  eddy  viscosity  in  incompressible 

> Sc  t 

flow  (p.^)  can  be  obtained  by  means  of  the  following  expression  which  has 
been  found  to  correlate  some  experimental  data  (Reference  2). 


p/  =  pK(26')(u  -  u  )  =  constant  ov  erase  ct  ion  (30) 

t  e  o 

where  K  is  on  the  order  of  0.02  to  0.005  in  axially  symmetric  and  two- 

5|< 

dimensional  flows,  and  6  is  the  semi-thickness  of  the  incompressible  wake. 
(It  should  be  noted  that  in  the  literature  the  quantity  p.^/  p  is  often  described 
by  the  terminology  "eddy-viscosity11  ,  rather  than  p.^  alone).  It  has  been 


welL  established  that  p/‘  >>  p.  in  regions  away  from  laminar  sub-layers. 


In 


* 


fact  p.^  is  on  the  order  of  hundreds  of  times  larger  than  p. '  in  wakes. 


Likewise  in  compressible  flow  it  is  expected  that  pi^»  pi. 

In  view  of  the  paucity  of  information  concerning  \ ^  and  (D_)  »  it 

has  been  considered  reasonable  to  assume  P  =  p.  c  /\  =P(  laminar) 

t  *t  p/  t 

Le  =  pD..  cpA  =  Le  ,  so  that  only  p.  will  require  explicit  evaluation, 
t  1J  t  t 
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Therefore,  in  a  rough  approximation,  if  turbulent  density 
fluctuations  are  negLected,  if  the  Laminar  transports  are  neglected  in 
comparison  with  the  turbulent  values,  and  if  it  is  assumed  that  P  =  P  (Lami¬ 
nar),  and  Le^  =  Le  (laminar),  the  Laminar  anaLysis  previously  given  can  be 
carried  over  directly  to  the  turbulent  case  with  the  provision  that  p  must 
replace  p  everywhere  except  in  p  ,  which  can  remain  intact  since  it  is  only 
a  reference  constant. 

It  is  observed  that  in  the  integral  method  analysis  p  appears  onLy  as 

p  in  the  transformed  variable  s  (Eq.  7)  and  in  the  chemical  rate  term  in 
o 

(10b).  It  has  been  pointed  out  that  the  chemical  rate  term  will  be  negligible 
with  respect  to  the  diffusion  term  in  (10b)  for  laminar  flow  at  the  altitudes 
of  immediate  interest  here.  Since  it  is  well  known  that  p^i^p*  and  there¬ 
fore  that  (perhaps  by  a  factor  of  severaL  hundred),  the  chemical 

processes  will  be  even  less  important  with  respect  to  diffusion  under 
turbulent  conditions.  This  is  physically  reasonable  since  mass  diffusion 
would  be  expected  to  be  enhanced  considerably  by  the  effect  of  turbuLence. 
Likewise  the  entire  streamwise  scale  of  the  wake  will  be  reduced 

in  roughly  proportion  to  l/p  .  This  means  a  streamwise  scale  reduction 

to 

by  a  factor  of  several  hundred.  That  is,  the  veLocity  parameter  u  wilL 

o 

approach  u  giving  a  uniform  profile  in  a  much  shorter  distance  in  turbuLent 
e 

flow  than  in  laminar  flow. 
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A  somewhat  rational  method  for  estimating  p.  p.  *  can  be  obtained 
on  the  basis  of  work  by  Mager  (Reference  16)  concerning  turbulent  compres¬ 
sible  boundary  layers  over  surfaces,  and  its  extension  to  the  case  of 
axially- symmetric  wakes  can  be  obtained  on  the  basis  of  work  by  Ting  and 
Libby  (Reference  7).  This  work  consists  of  stating  a  set  of  conditions  which 
are  sufficient  to  transform  the  turbulent  compressible  equations  of  momentum 
and  overall-continuity  to  their  incompressible  form  .  The  transformation 
of  Mager  has  been  utilized  with  some  success  to  correlate  experimental 
results  in  flow  over  bodies.  Its  primary  conditions  are  that  the  Howarth- 
Dorodnitsyn  density-transformation  of  the  normal  coordinate  can  be  intro¬ 
duced,  that  the  stream -function  remains  invariant  under  the  transformation, 
and  that  the  turbulent  shear  within  a  mass-element  remains  invariant  under 
the  transformation.  The  latter  condition,  which  is  by  no  means  obvious, 
actually  stems  from  a  physical  interpretation  of  one  of  the  mathematically 
required  conditions  for  reduction  of  the  equations.  Another  aspect  of  the 
transformation  requires  the  somewhat  arbitrary  selection  of  values  for 
p  and  p.  ‘  in  the  incompressible  flow. 

The  same  conditions  are  required  for  the  axially-symmetric  wake 
except  that  now  the  invariance  of  the  moment  (about  the  axis)  of  the  turbulent 
shear  within  a  mass  element  is  required  for  a  reduction  of  the  equations. 

The  validity  of  the  transformation  for  the  axially  symmetric  wake  still  awaits 
experimental  check.  In  this  case  Ting  and  Libby  have  suggested  that  the 


I*— 
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p*  and  |i  *  be  evaluated  at  conditions  at  the  outer  edge  of  the  wake. 

Neither  Mager  nor  Ting  and  Libby  derived  the  relations  required 

to  reduce  the  energy  equation  and  the  species-conservation  equations  to 

incompressible  form.  These  wouLd  yield  relations  for  X  /  \  *  and 

(D..)  /  (D..)*.  It  would  be  interesting  to  examine  the  consistency  of  such 

XJ  t  ij  t 

relations  and  the  assumptions  P  =  P  and  Le  =  Le.  This  will  be  Left  to 

l  L 

future  work.  An  assumption  in  the  turbuLent  energy  equation  that  the 


turbulent  stagnation-entha Lpy  transport  v  'H'is  proportional  to  the  momen¬ 
tum  transport  u^1  (primes  denote  turbulent  fluctuations  and  bars  mean 
quantities)  Leads  to  the  isobaric  Crocco  reLation  between  mean  variables. 
That  is,  H  =  A  +  Bu,  where  A.  and  B  are  constants.  The  condition 


v'H1  -  Bu'v1  is  satisfied  when  P^  *  1  and  H  =  A .  +  Bu.  The  special 
case  for  which  B  =  0  corresponds  to  isoenergetic  mean  flow. 

The  following  reLation  was  derived  by  Ting  and  Libby  from  the 
condition  of  invariance  of  the  turbulent  shear-moment: 

H-t  /  Kt*  =(1/py2)J0  pd(y2) 

It  may  be  assumed  from  (30)  that 

i a  *=  p  K(26  )(u  -  u  ) 

t  e  m  e  o 

In  our  integraL  method  anaLysis  wiLL  require  evaluation  onLy  at 
the  symmetry  axis,  that  is, 

at  y  =0,  n  =  u  =  u  *  =  p  K  (26  )(u  -  u  )  (31c) 

t  to  t  e  m  e  o 

In  some  cases  a  constant  mean  vaLue  of  p.  *  may  be  used  with 
sufficient  accuracy. 


(31a) 


(31b) 
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III  CALCULATIONS 

CaLcuLations  have  been  made  of  the  far-fieLd  Laminar  wake 
properties  at  the  foLLowing  trajectory  points  of  a  bLunt  body: 


Code 

Altitude  (ft) 

VeLocity  (fps) 

1 

150, 000 

18, 900 

2 

200, 000 

20, 200 

3 

250, 000 

20, 400 

The  corresponding  initiaL  conditions  at  the  axis  of  the  wake  and 

other  pertinent  quantities  are  shown  in  TabLe  I.  It  is  recaLLed  that  the 

initiaL  conditions  are  obtained  by  assuming  that  chemicaL  equilibrium  is 

achieved  reasonably  cLose  behind  the  bow  shock  and  that  the  fLow  then 

expands  in  chemicaLLy  frozen  fashion  but  with  vibrationaL  equilibrium  untiL 

ambient  pressure  is  achieved.  This  expansion  is  isentropic.  Insofar  as 

eLectron  production  is  concerned,  the  same  procedure  is  empLoyed  for  the 

NO*  reaction,  utilizing  (28)  to  give  the  equilibrium  vaLues,  and  then  simpLy 

muLtipLying  the  frozen  mass  fraction  ajqQ+  the  initiaL  wake-density 

p  and  appropriate  conversion  factors  to  get  the  initiaL  number -density  of 
oc 

NO*  and  therefore  of  e”. 

Consideration  of  the  relations  (10b)  and  (  24)  for  species  concen¬ 


trations  of  O  and  N  in  the  wake  indicate  immediately  that  the  chemicaL  rate 
production  term  is  negLigibLe  with  respect  to  the  diffusion  term  for  the  fLight 
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conditions  Listed  above.  However,  this  is  not  the  case  for  the  NO  and  e " 
reaction.  Initially,  the  ionization  recombination  rate  term  dominates  and  is 
shown  to  cause  a  rapid  decrease  in  the  number-density  of  electrons.  For 
example,  for  all  flight  conditions,  the  number-density  of  electrons 
(particles/cc)  decreases  by  at  least  a  factor  of  one  hundred  within  a  distance 
of  ten  feet.  This  strongly  indicates  that  the  electron  species  is  not  frozen 
in  the  anterior  region  of  the  wake,  and  with  the  extremely  high  rate  of 
recombination  one  should  suspect  that  the  electron  species  in  the  wake  is 
close  to  chemical  equilibrium. 

It  should  be  noted,  that  for  chemical  equilibrium  of  the  NO^  or  e 
species  the  electron  density  is  properly  obtained  only  by  setting  the  bracket 
in  equation  (28a)  equal  to  zero. 

Figures  8  through  14  present  the  results  of  specific  calculations. 

In  particular,  the  variations  of  quantities  along  the  axis  are  given,  where 
the  independent  parameter  is  in  the  physical  length  (feet).  For  each  tra¬ 


jectory  point  calculations  were  made  under  two  assumptions 


P  =  Le  =  1 
o  o 


and 


P  =  0.71,  Le  =  1.45 
o  o 


for  the  isoenergetic  initial  conditions  H  =1, 

oc 


and  for  the  200,  000  feet,  20,  200  fps  case  calculations  were  made  with 


P  =  Le  =  1 
o  o 


for  an  initial  defect  in  H  ,  that  is,  H  =0.9. 

o  oc 


The  stagnation  enthalpy  His  shown  in  Figure  8a  for  the  isoenergetic 

initial  condition.  It  is  interesting  to  note  that  the  maximum  defect  in  H  for 

o 


the  non -isoenergetic  wake 


P  =  0.71,  Le  =  1.45,  H  =  1 
o  o  oc 


is  about 
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7  to  9  percent  and  agrees  cLoseLy  with  the  prediction  of  Appendix  B.  Shown 

in  Figure  8b  is  the  variation  of  H  at  200,000  feet  as  a  resuLt  of  an  initial. 

o 

defect  in  H  ,  that  is,  H  =0.9. 
o  oc 

The  veLocity  decay  aLong  the  axis  is  presented  in  Figures  9a  and  9b 

for  the  isoenergetic  and  non-isoenergetic  initial,  condition  respectively. 

Figures  10a,  11a,  12a  respectively  illustrate  the  variations  of  mass  fraction 

of  air  (  a  ),  temperature  (  T  )  and  density  (  p  )  along  the  axis  for  the 
ao  o  o 

initial  isoenergetic  condition  H  =  1.0.  The  Prandt L-Lewis  number  effect 
fe  oc 

is  clearly  illustrated  by  these  curves.  Of  main  importance,  is  the  substan¬ 
tial  increase  in  the  maximum  value  of  the  temperature  for  the 

case  over  the  isoenergetic  wake;  especially  if  one  is  interested 
in  the  radiation  properties  of  the  wake.  Analogously,  the  above  variations 


Le  =1.45 
o 


P  =0.71, 
o 


for  an  initial  defect  in  H  ,  that  is,  H  =0.9  are  present  in  Figures  10b,  lib, 

o  oc 


12b  for  the  200,  000  feet,  Jj3o=  L. eQ=  lj< 


case. 


The  mass  fraction  of  NO  along  the  axis  as  calculated  from  (10b)  and 

(28)  is  shown  in  Figures  13a,  b,  c  for  the  isoenergetic  initial  conditions 

H  =1  and  Figure  13d  presents  the  variation,  for  the  200,000  feet,  P  =  Le  =1 
oc  o  o 

case,  with  an  initial  defect  in  H^.  The  corresponding  number -density  of 

e 

electrons  (nQ)  are  shown  in  Figures  14a,  b,  c,  d.  F rom  these  figures,  the 

effect  of  an  initial  defect  in  H  can  be  determined.  Figure  14b2  shows  that 

o 

the  initial  rate  of  decrease  of  the  number-density  is  larger  for  the  non-iso- 
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energetic  initial,  conditions. 

FinaLLy,  the  physical.  scaLe  x  of  the  streamwise  coordinate  which 
is  reLated  to  s  by  (7)  is  shown  in  Figures  15a  and  15b.,  Notice^that  the 
effect  of  aLtitude  is  to  increase  the  physical.  Length  at  the  Lower  aLtitudes 
and  decrease  it  at  higher  ones.  This  is  due  to  the  dependence  of  the  physicaL 
Length  (x)  on  the  free  stream  ReynoLds  number.  It  can  aLso  be  seen  that 


PrandtL-Lewis  number  effect 


P  =  0.71,  Le  =  1.45 
o  o 


causes  a  slight 


P  =Le  =1 
o  o 


decrease  in  the  physicaL  Length  over  the  isoenergetic  conditions 
This  is  due  to  the  inverse  dependence  of  the  physicaL  Length  (x)  on  temperature 
For  turbulent  flow,  the  foLLowing  expression  may  be  written 
from  (31c): 


'to  


Pe  ue 


Ho 


1/2 

2K  (  6m  ]  .  ,  -  . 

(l’Uo> 


(32) 


In  our  present  calculations,  this  quantity  varies  from  the  order 
3  5 

of  L0  at  250,000  ft.  to  10  at  L50,000  ft.  This  gives  an  idea  of  the 
substantial  relative  decrease  in  the  length  of  the  wake  in  turbuLent  flow, 
based  on  the  assumptions  previously  stated.  Under  these  assumptions, 
at  250,000  ft.  the  turbuLent  wake  is  on  the  order  of  L000  times  shorter  than 
the  laminar  wake;  while  at  L50,000  ft.  the  turbuLent  wake  is  about  100,000 


times  shorter  than  the  Laminar  wake. 
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A 

B 

C 


c 

P 

D 


D. 

1 

Dij 


E 


E 


h 

P 

H 

I. 

1 

k 


k 

r 

K 


Le 


NOTATIONS 


constant 

constant  (aLso  see  Eq.  Bl) 

constant  (Eq.  29) 

specific  heat  ;  constant  pressure 

diameter 

dissociation  energy  per  mole  associated  with  atomic  species 
binary  diffusion  coefficient  for  the  interdiffusion  of  species  i  and  j 
partition  function 

electronic  partition  function  of  atomic  species 
electronic  partition  function  of  molecular  species 
static  enthalpy 
Plank's  constant 
total  enthalpy 

ionization  energy  per  mole  associated  with  ionized  species 

thermal  conductivity  coefficient 

forward  reaction  rate  parameter 

reverse  reaction  rate  parameter 

defined  by  eq.  30 

pD  •  -c 

Lewis  number  Le  -  P 


k 
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m 

m. 


1 

M 

1 

M 


n 

i 

N 

N2 

O 

o2 

P 

P 

R 


s 


T 


V 


T 

i 


u 

u 


>  v 

oo 


V 

V 
i 

w 


transformed  normal  coordinate 
particLe  mass 

moLecuLar  mass  mass/moLe 

reduced  mass  (Eq.  24e) 

free  stream  (flight)  Mach  number 
,  th 

particLe/cc  of  i  species 

transformed  normal  coordinate  except  where  otherwise  stated 

moLecuLar  nitrogen 

atomic  oxygen 

moLecuLar  oxygen 

pressure 

p  cn 

PrandtL  number  P  -  _ tL 

k 


gas  constant  per  mole 

transformed  streamwise  coordinate  (aLso  see  Eq.  23a) 
temperature  °R 

characteristic  vibrationaL  temperature 
streamwise  and  normal  veLocity  component 
free  stream  (flight)  veLocity 
LocaL  mean  veLocity 

t  h 

diffusion  veLocity  of  i  species 

th 

net  production  of  i  species  (L/sec) 


i 
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x,  y 

streamwise  and  normal  coordinate 

X 

cataLyst 

a  . 

1 

t  h 

mass  fraction  of  i  species 

6 

boundary  Layer  thickness 

6 

m 

transformed  boundary  Layer  thickness 

6* 

semi-thickness  of  incompres sibLe  wake 

0 

momentum  thickness  0"^opu(l-  u)  ydy 

0 

E 

energy  thickness  0  =  J  pu  (1  -  H)  ydy 

^  ° 

9 

i 

th 

i  species  thickness  0*U=JO  p  u  (a.^  -  a  ydy 

P 

absoLute  viscosity 

P 

density 

Pdi 

characteristic  density  for  moLecuLar  species  (Eq.  24f) 

A 

defined  (Eq,  2  5) 

SUBSCRIPTS 


c 

initial  conditions  for  wake  caLcuLations 

e 

conditions  outside  viscous  layer 

e~ 

refers  to  eLectron 

i 

£  .  th 

refers  to  i  species 

N 

refers  to  atomic  nitrogen 

n2 

refers  to  molecular  nitrogen 

NO+ 

refers  to  ionized  nitric  oxide 
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o  vaLues  evaluated  along  axis 

O  refers  to  atomic  oxygen 

Oz  refers  to  molecular  oxygen 

s  denotes  total  differentiation  with  respect  to  indicated  variable 

st  stagnation  values  behind  normal  shock 

t  turbulent 

x,  y  denotes  partial  differentiation  with  respect  to  indicated  variable 
oo  undisturbed  flight  conditions 


SUPERSCRIPTS 


denotes  non-dimensional  quantities  unless  otherwise  specified 


u 

u  = —  , 


p  —  H  -_IJL  h 

Te'  H  *  H  '  h  =  T( 


T  = 


* 


incompressible 
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APPENDIX  A 

The  equations  governing  the  flow  of  multicomponent  gas  mixtures 
are  discussed  in  considerable  detail  by  HirschfeLder,  Curtiss  and  Bird 
(Reference  14;  see  particularly  Chapters  7,  8  and  11).  In  this  appendix 
we  shall  pay  particular  attention  to  the  diffusion  velocities  Vi  of  the 
several  species.  The  diffusion  velocity  of  a  given  species  represents  the 
mass -averaged  velocity  of  the  species  measured  with  respect  to  the  mass- 
averaged  velocity  of  the  mixture.  Here,  the  mixture  velocity  is  represented 
by  the  components  u  and  v. 

The  boundary -Layer  equation  of  mass  conservation  of  species  in 
the  wake  is  given  in  terms  of  V.  as  follows: 

( puy  a  •)  +  ( pvya .)  =  (ypa.V.)  +  pyw.  (A.  1) 

1  x  1  y  l  1  y  l 

At  the  symmetry  axis,  y  =  0,  this  equation  takes  the  form 

ua  =  2a  .V.  +  pw.  (A .  2) 

Contributions  to  the  diffusion  velocity  are  made  by  several  pheno¬ 
mena.  Only  one  cLass,  gradients  of  the  concentration  of  species,  are  ex¬ 
pected  to  be  significant  in  the  present  instance.  Their  effects  are  usually 
termed  "mass  diffusion".  The  other  processes,  such  as  thermal  diffusion, 
pressure  diffusion,  externaL  forces,  and  the  Dufour  effect,  are  neglected 
here.  The  diffusion  velocities  in  a  gas  with  v  components  can  be  expressed 
in  terms  of  the  binary  -  diffusion  coefficients  D^.  (not  to  be  confused  with  the 
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more  artificial  "multicomponent  diffusion  coefficients")  by  the  following 
set  of  relations  (Reference  14,  p.  516)  : 


(l/n)  “  Vi>  *  Vin(n{  /  n) 


(A  .  3a) 


j  =1 

j/i 


^  14  nr^  V*  *  0 

i 

where  m  is  a  number  density  (particles/cc) ,  n  =  ^  m  ,  m  rm  =  p  a. 

i 

m.  is  the  particle  mass,  and  D..  =  D... 
i  F  ij  Ji 

For  a  binary  mixture  (A.  3a)  yields  one  relation  and  (A.  3b) 

another.  From  these  we  obtain 

V.  =  -D.  ,  yin  a  .  i  =  1,2 
i  i 2  7  i 


(A.  3b) 


(A.  4) 


thus  indicating  the  significance  of  the  binary  coefficients  D  — •  Theoretical 
and  experimental  data  for  estimating  the  transport  properties  D_  for  a 
variety  of  pairs  of  species  are  given  in  Reference  14. 

For  other  than  binary  mixtures  the  expressions  for  the  V\  are  more 
complex.  In  particular  we  shall  be  interested  here  in  a  quadruple  mixture 
of  O,  N,  Oz  and  N2  .  In  generaL  each  V\  is  a  Linear  function  of  the  concen¬ 
tration  gradients  of  the  several  species,  the  coefficients  being  functions  of 
the  D  and  the  species  concentrations.  However,  when  the  particle  weights 

ij 

of  the  species  are  not  extremely  different  from  each  other,  the  D..  will  not 

ij 


GENERAL  APPLIED  SCIENCE  LABORATORIES,  INC. 


PAGF  46 


differ  substantially  in  magnitude  from  each  other.  Species  such  as 
O,  N,  02  ,  N2  and  NO  can  be  considered  suitable  for  approximations 
involving  average  particle -weights  and  average  diffusion  coefficients. 

An  approximation  along  these  lines  will  now  be  suggested  for  use  here. 

The  quantity  V  in  (m  /  n)  is  re-expressed  in  terms  of  the  mass 
fractions  q.  as  follows: 


Vln(ni/„»  “(„)(  a) 

V  QJ  Yai\Vaj 

Z_  mj  Vai'  mj 

j/i  _ 

(A.  5a) 


For  a  binary  mixture,  the  term  enclosed  in  square  brackets  in 

(A.  5a)  is  equal  to  l/m..  Likewise,  for  an  average  value  of  the  particle 

J 

weights  m.,  but  for  any  number  of  components,  the  term  in  brackets 
J 

equals  1  /  (m.)  ;  in  this  case, 

J  av 

Vin(n.  /n)  -(l/m.)  (p/n)(  VQ  /a.)  (A. 5b) 

i  J  av  1  1 

The  left  side  of  (A.  3a)  can  be  rewritten  as 
v 


(p/n)  ^  Uj/mjDij)(vj  -  v.) 
j  /  i 


(A.  5c) 


For  average  values  (m.)  and  (D..)  this  becomes 

J  av  ij  av 


( p/n)(  l/m  .D . .)  (a  V.  -  V  +  aV)  ^(p/n)(l/mD  )  (-V)(A.5d) 

J  ij  av  1  1  i  ii  j  ij  av  i 
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Equation  of  (A  .  5b)  and  (A .  5d)  yields  an  approximation  to  (A  .  3a) 
for  a  muLticomponent  mixture  with  particle  weights  of  the  same  order, 
name  ly, 

V.  -  -(D  .)  Va./a  (A.5e) 

1  ij  av  i  i 

This  is  the  expression  used  in  (4)  with  Le  3  p(D  .)  c  /\, 

ij  av  p 

P  =  pc  / \,  c  =\  a  . c  .  The  V.  obtained  from  this  relation  automatically 
P  P  i  Pi  i 

satisfy  (A  .  3b) . 

With  (A.5e)  the  heat  flux  which  appears  in  the  energy  equation  (3) 
takes  the  following  form 


q  =  -  WT 


+  p  a- 


.  h.  V. 

Ill 


=  -(p/P)Vh+  (ji/P)(l  -  Le)^T  h.Va. 


(A.  6) 


Although  there  would  be  no  difficulty  in  principle  in  obtaining 

solutions  with  more  exact  representations  of  the  multicomponent  V  ,  the 

i 

computational  program  would  be  substantially  longer.  Therefore  it  was 
considered  to  be  sufficient,  at  least  for  initial,  calculations,  to  use  (A.5e) 
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APPENDIX  B 

Estimate  of  PrandtL  number  -  Lewis  number  effects  on  the  stagnation 

enthaLpy  parameter. 

Consider  (19c)  in  the  following  form,  with  H  =  H  and  u  =  u 


dH/du  =  (1  -  H)  /(I  -  u)  +  B 


(B.l) 


where 


B  = 


1  -P, 


1-H  u( 
1  -u  "  H 


u 


e  _ 


(Le  -  1)  ^Lh.  (  a .  -  a  ) 

o  i  io  le  10 


H  (1-u) 

e 


The  summation  of  static  enthaLpy  in  B  may  be  treated  as  foLLows 


h.  (a .  -  a .  )  = 

/  io  le  io 


Zh.  a-  -h  SY  h-  a.  -  (  H  -I*-) 
io  le  °  /  io  le  o  2 


(B.2) 


In  the  air-wakes  under  consideration  here,  it  is  assumed  that  the 

atomic  species  in  the  outer  fLow  are  negLigibLe.  Hence  the  a  terms 

ie 

represent  onLy  moLecuLar  species.  According  to  the  conventions  of 
statistical  thermodynamics,  the  dissociation  energies  do  not  appear  in  the 
expressions  for  the  enthaLpies  of  the  moLecuLar  species.  Hence,  for  the 


moLecuLar  species  treated  as  perfect  gases 

fT 

h.  =  J  c  (T)  dT 

1  °  Pi 


(B .  3) 


where  c  takes  into  account  the  trans LationaL  and  vibrational,  degrees  of 
Pi 

Z  ,  7  O 

freedom.  For  02  and  N2  ,  Cp  has  a  vaLue  on  the  order  of  8000  ft  /  sec  R 
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at  elevated  temperatures.  In  this  estimate  we  need  not  distinguish 

between  Oz  and  N2  in  view  of  their  similar  molecular  weights.  Hence 

we  may  consider  a.  =  1  and  h. ~ (8000  ft2  /sec2  °R)  T,  and  (B.2)  becomes 

ie  i 

(8000  ft2/sec2  °R)T  -  (H  -  u2  /  2)  (B.3) 

o  o  o  ' 

At  hypersonic  speeds,  u*  /  2,  and  division  of  (B.  3)  by 

H  yieLds 
e 

(8000  ft Z/  sec 2  °R)T  /(u  2  /  2)  ~  (H  -  u2  )  (B.4) 

o  e  ' 

F or  u^  on  the  order  of  20,  000  ft/ sec  (B.  4)  becomes 

0.04(  T/1000°R)  -  u2  (B.5) 

In  the  far-field  wake  H  is  expected  to  have  values  on  the  order 

of  unity,  u  varies  between  approximately  0.5  and  1.0,  and  temperatures 

on  the  order  of  several  thousand  °R  may  be  expected.  Hence  the  factor 

in  (B.5)  will  be  a  negative  fraction,  and  the  corresponding  coefficient  of 

(Le  -  1)/P  in  the  definition  of  B  wiLL  be  a  negative  quantity  on  the  order 

of  unity.  This  estimate  probabLy  is  valid  even  as  u~ >1,  since  the  numerator 

of  the  coefficient  as  represented  by  (B.5)  becomes  very  small  in  this  regime. 

The  coefficient  of  (1  -  P  )  /P  in  the  definition  of  B  may  be  assumed 

o  o 

with  the  use  of  (18)  under  hypersonic  conditions  as  follows: 

( 1  -H)/(  1  -u)  -  (u  2/H  )  u  *  (1-H  )/(l-u  )  -  2u  - 

e  e  c  c 


-  2u 


(  B.  6) 
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Therefore 

(B.7) 

where  l|j  is  a  positive  quantity  on  the  order  of  unity.  Since,  for  air 

1-P  ^  Le  -  1,  it  is  seen  that  for  the  present  conditions  of  interest,  B 

o  o 

is  a  negative  quantity  whose  magnitude  is  on  the  order  of  (1  -  P  )  /  PQ. 

For  a  constant  average  value  of  B,  Eq.  (B.  1)  may  be  integrated 

to  yield 


B  -  - 


(1-Po>  „  .  ,  l^o-D 

-  (2u)  + - i|j 


‘-H  ‘‘He  .  ,  /I '“A 

: -  *  - -  +  /B/  Jnl  - - 1  (B.l 

1  - u  1-u  \ 1 - u  / 

c 

where  the  integration  constant  has  been  evaluated  at  the  initial  station 
where  values  are  designated  by  the  subscript  c. 

Thus  it  is  seen  that  the  Prandtl  number  -  Lewis  number  term  has 


the  effect  of  delaying  the  approach  of  H  to  unity  as  u~->l. 


TABLE  I 
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150, 000' 

200, 000' 

250, 000' 

FLIGHT 

CONDITION 

18,  900 '/ sec 

20, 200 '/ sec 

20, 400 '/sec 

Pco 

#/sq  ft 

2.940 

0.460 

0.040 

P°0 

slug/ft3 

3.39  x  10"6 

5.88  x 

io“7 

,  -8 

6.  55  x  10 

T 

OO 

^  CO 

°R 

#  sec/ft2 

504 

-8 

35.8  x  10 

454 

-8 

33.8  x  10 

354 

-8 

27.  1  x  10 

M 

oo 

He=  hst 

17.2 

19.3 

22.3 

ft 2  /sec  2 

8 

1.81  x  10 

2.07  x 

8 

10 

8 

2. 10  x  10 

pst 

#/ sq  ft 

1210.0 

240.0 

27.  30 

Pst 

slug/ft 3 

-5 

4.  67  x  10 

9 . 2  6  x 

-6 

10 

-6 

1. 11  x  10 

T 

st 

ast 

°R 

11,000 

10, 400 

9,  700 

mass  fraction 

0.3  65 

0.445 

0.475 

a° 

st 

mass  fraction 

0.232 

0.232 

0.232 

a  N 

St 

mass  fraction 

0.133 

0.213 

0.245 

NO+ 

Qst 

mass  fraction 

,  -3 

0.563  x  10 

0.566  x 

-3 

:  10 

-3 

0.485  x  10 

NO+ 

nst 

Pco 

particles 

cc 

slug/ft3 

2. 73  x  1014 
-7 

7.42  x  10 

5.45  x 

1 . 35  x 

13 

10 

-7 

10 

12 

5 . 60  x  10 

1.49  x  10"8 

T 

CO 

°R 

1672.0 

1352.  0 

1125.  0 

O  O 
+ 

particles 

CC 

4.44  x  1012 

8.13  x 

11 

10 

10 

6. 89  x  10 

H 

CO 

1.0 

1.0 

0.9 

,  1.0 

\ 

u 

CO 

0.  692 

0 . 65  6 

0.574 

0.  634 

e 

CO 

198.  6 
C  =  — pp — 

0.0065 

0.394 

0.0071 

0.437 

0.0083 

0.0074 

0.564 

oo 

\ 

i 

i 

XNO+ 

electron  volt 

9.258 

9.258 

9.258 

dno 

electron  volt 

6.42 

6.42 

6.42 

Do2 

electron  volt 

5. 1155 

5.1155 

5. 1155 

l 

Dn2 

electron  volt 

I  9.756 

9.  756 

| 

f 

9.  756 

j 

1 
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Figure  1.  Schematic  Diagram  of  Far-Field  Wake 
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Figure  2c-  Initial  Mass  Fraction  of  Atoms  for  Frozen  Flow 
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Figure  3.  Quadratic  Velocity  Profile  (Eq.  13a) 
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Plot  of 
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versus  uQ 


Figure  4. 


(Eq.  16) 


Figure  5a.  Velocity  parameter  versus  streamwise  coordinate  s  (Eq.  17) 


Figure  5b.  Velocity  parameter  versus  streamwise  coordinate  s  (Eq.  17) 
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Figure  6.  Rough  Estimate  of  Effects  of  Prandtl  Number  and 

Lewis  Number  Deviations  from  Unity  on  the  Stagnation 
Enthalpy-Velocity  Relation. 
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Figure  7a.  Concentration  parameter  versus  velocity  for  diffusion-controlled 

species  with  Le  /P  =1  and  2. 
r  o  o 


Figure  7b.  Concentration  parameter  versus  streamwise  coordinate  for  diffusion-controlled  species 
with  Leo/PQ  =  1  and  2.04. 


GENERAL  APPLIED  SCIENCE  LABORATORIES,  INC. 


PAGE 


65 


o 


Figure  8b.  Stagnation  enthalpy  variation  aLong  the  axis. 


Figure  9b.  Velocity  Decay  Along  the  Axis 
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Figure  10b.  Mass  fraction  of  air,  variation  along  the  axis. 
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Figure  lib.  Temperature  of  air,  variation  along  the  axis 
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Figure  12b.  Density  of  Air,  Variation  Along  the  Axis 
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Figure  13c. 
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Figure  14d.  Variation  of  Electron  Density  Along  the  Axis 
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Figure  16.  Variation  of  Electrons  in  the  Inviscid  Flow  About  a  Blunt  Body.  FIGURE  16  (Revised) 
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SUPPLEMENT  TO  GASL  TECHNICAL  REPORT  NO.  180 


HYPERSONIC  AXISYMMETRIC  WAKES  INCLUDING  EFFECTS  OF 

RATE  CHEMISTRY 
By  M.  H.  Bloom  and  M.  H.  Steiger 

Retransformation  to  physical  coordinates  by  use  of  the  approximation 
1/p  =  h. 

The  inverse  transformation  is  performed  by  means  of  the  following 

equations  see  (7)  : 

L  J  /N 

y2/2=  62  J0  ( l/p )  NdN 


L  U/P' 


6  /2  =  6m'/0  (1/P)NdN 


We  shaLL  assume 


Now 


and  from  (13a) 


where 


l/p  =  h  where  h  =  h/h  ,  p  =  p/p 

e  6 


h  =  H  -  u2  /2 


H  =H  +  (H  -  H  )  f 
o  '  e  o 

u  =  u^  +  (u  -  u  1  f 
o  '  e  o' 

2  3  4 

f  =f (N)  =  6N*-  -  8N  +  3N 


Thus 


fi  =  1  -  (l-H0)(l-f)  +  (ue2/2hj  -j 


2(l-%)-(l-H0) 


(1-  f)  -  (1-uq)2  (1-f)2 
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Now  define 
N 


=  fQ  (1  -  f )  NdN  =  -y-  -  |  N4  +  |  N5  -  4  N6 


i 


N 


h  =Jo  ll-£»  NdN 

=  y  N2  -  3N4  +  —  N5  +  5N6+8N8+  2:N10-  ?rN7+f  n8-  4 r  N? 
&  d  ID  7  2  3 


he  =L  ('-£)  NdN  =  jo 


2e 


( l-f )  NdN  = 


11 

Z10 


Thus 


(6/6m)2=1  -2Ile(£-H0)^(u£/2he)|[a(l-ilo)-(l-H0)]  Iie-(l-S0)2Iie 


In  the  calculations  of  GASL  TR-180  an  initial  value  of 


6  =  0. 5  ft 

me 


was  used  for  all  flight  conditions. 


F or  the  special  case  of  Hq  =  1 ,  we  derive 


(6/6m)2  =1  +  (u2e  /2he  )(2A)(l-uo) 


1  -  (1 1/4Z)(1  -uo)j 


where  6  2 
m 


21.0 

T r 

o”  10 


(i-ujU  +-U-  aQ) 
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giving 


and 

6 


r  f  ..2'  x”i 

i 

8.46  fne\ 

c  l2irJ 

2 

L  \  e  - 

1 


U 

42 


(i-u  )  + 
o 


i 

1/2 

1 

1/2. 

i  (i)  d-uj 
2hg  K5  '  °  J 

1  +  1.10  u 

o 

This  reLation  is  pLotted  in  the  attached  figure  for  the  typical.  vaLue 

(2/5)(u2  /2h  )  =32.  Since  h  =c  T„  at  low  T  ,  u2  /2h  =  Ye_1  m2  /5, 

e  e  e  pe  e  e  e '  e  - £ -  e 

and  this  corresponds  to  Me=20. 

It  is  noted  that  6  is  approximately  constant  for  uQ  <  0.95.  This 
constitutes  the  major  portion  of  the  wake  dealt  with  here  (see  Figures  9a,  b). 
Moreover  the  factor  which  nondimensionalizes  6  has  the  approximate 


value 


=  0.11  Me  (ft) 


.'.  6  =a(0. 11  M  )  for  u  <0.95 
v  e  o 

where  MaM  is  a  constant  factor  estimated  from  the  aforementioned  plot. 

In  this  case,  a  =  0.8. 

If  0.9  <  H  <  l  the  effect  on  6  of  deviations  from  unity  will  be  small, 
o 

Likewise  the  changes  in  "a"  for  Mach  numbers  M^  on  the  order  of  20  will  not 
be  strongly  influenced  by  deviations  in  Me  from  20. 


M.  H.  Bloom 
January,  1961 


.108  Me  (ft) 
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25 


Eq.  (17c)  f(uQ)  i  (5/l2)(l/l-u0)  -(iy6>in 


(10  +  11  uQ)/(l  -  uQ) 


/  m.  \  i/2  “T.  /  T  r'  /  v' 

Eq.  (24f)  Pdi  *\J^~ )  (,,mikT)  (Ti  /T)(‘  -  * 


Eq.  (25)  h  =  RT 


5  faO  .  °n\  aP2  ( 7  /  \  QN2  ( 7  a  ^ 

2  U  +MJ  +  Mo,  v  °2/+  mn,V  Y 


(  DQ2  \ 

+  QoV  2mo/  +  QN  \2mN 


Eq.  (28b)  kj  .  /  k  , 
dNO+/  rNO+ 


MNO+  Me  | 

mnmo 


3/2 


10.42 


6111 
1  -  e  '  “T“ 


."^NO  "  °NO^/ 
e 


kT 


„  -41  IK  -589K 
5  +  3e  +  e 
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APPENDIX  C  TO  GASL  TECHNICAL  REPORT  NO.  180 

HYPERSONIC  AXISYMMETRIC  WAKES  INCLUDING  EFFECTS  OF 

RATE  CHEMISTRY 
By  M.  H.  Bloom  and  M.  H.  Steiger 

The  attached  figures  present  the  results  of  the  analysis  described  in 
GASL  Technical  Report  No.  249,  Improved  Hypersonic  Laminar  Wake  Cal¬ 
culations  Including  Rate  Chemistry,  by  Martin  H.  Steiger. 

Axial  variations  of  temperature,  density,  electron  density  and  velocity 
deficiency  are  given  for  the  flight  conditions  listed  in  Table  I. 

Free  stream  pressure,  temperature  and  density  have  been  revised  in 
accordance  with  the  Air  Force  Cambridge  Research  Center  Atmospheric 
Model  of  December,  1959. 
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